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Abstract 



The results from Suzaku observations of the galaxy cluster Abell 2256 are presented. 

^ ! This cluster is a prototypical and well-studied merging system, exhibiting substructures 

(-H . both in the X-ray surface brightness and in the radial velocity distribution of member galax- 

ies. There are main and sub components separating by 3!5 in the sky and by about 2000 

P ■ km s~^ in radial velocity peaks of member galaxies. In order to measure Doppler shifts of 

c/3 ■ iron K-shell lines from the two gas components by the Suzaku XIS, the energy scale of the 

' 

instrument was evaluated carefully and found to be calibrated well. A significant shift of 
the radial velocity of the sub component gas with respect to that of the main cluster was 
^ . detected. All three XIS sensors show the shift independently and consistently among the 

^ . three. The difference is found to be 1500 ±300 (statistical) ±300 (systematic) km s"^. The 

CN ■ X-ray determined absolute redshifts of and hence the difference between the main and sub 

^ ■ components are consistent with those of member galaxies in optical. The observation indi- 

cates robustly that the X-ray emitting gas is moving together with galaxies as a substructure 
within the cluster. These results along with other X-ray observations of gas bulk motions 

, , , in merging clusters are discussed. 

rS ■ 

. Key words: cosmology: large-scale structure — galaxies: clusters: individual 



(A2256) — intergalactic medium — X-rays: diffuse background 



1. Introduction 

Galaxy clusters are the largest and youngest gravitationally bound system among the hier- 
archical structures in the universe. Dynamical studies of cluster galaxies have revealed that some 
systems are still forming and unrelaxed. X-ray observation of the intracluster medium (ICM) pro- 
vided further evidence for mergers through spatial variations of the gas properties. Remarkably, sharp 
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X-ray images obtained by Chandra have revealed shocks (e.g. Markevitch et al. 2002) and density 
discontinuities (or "cold fronts"; e.g. Vikhlinin et al. 2001). These are interpreted as various stages 
of on-going or advanced mergers and suggest supersonic (for the shock) or transonic (cold front) gas 
motions. Cluster mergers involve a large amount of energies and hence influence numerous kinds of 
observations. In particular, possible effects of gas bulk motions on the X-ray mass estimates have 
been investigated extensively mostly based on numerical simulations (e.g. Evrard et al. 1996, Nagai 
et al. 2007). This is mainly because that the cluster mass distribution is one of the most powerful 
tools for the precise cosmology. Furthermore, cluster mergers heat the gas, develop gas turbulence, 
and accelerate particles, which in turn generate diffuse radio and X-ray halos. 

To understand physics of cluster mergers, gas dynamics in the system should be studied. The 
gas bulk motion can be measured most directly using the Doppler shift of X-ray line emission. These 
measurements are still challenging because of the limited energy resolutions of current X-ray instru- 
ments. Dupke, Bregman and their colleagues searched for bulk motions using ASCA in nearby bright 
clusters. They claimed detections of large velocity gradients, such as that consistent with a circular 
velocity of 4100t^?[!J] km s ^ (90% confidence) in the Perseus cluster (Dupke & Bregman 2001a) and 
that of 1600 ± 320 km s~^ in the Centaurus cluster' (Dupke & Bregman 2001b). These rotations 
imply a large amount of kinetic energy comparable to the ICM thermal one. Note that they used the 
ASCA instruments (GIS and SIS) which have gain accuracies of about 1% (or 3000 km s^^). Dupke 
& Bregman (2006) used also Chandra data and claimed a confirmation of the motion in the Centaurus 
cluster. These important results, however, have not yet been confirmed by other groups. For exam- 
ple, Ezawa et al. (2001) used the same GIS data of the Perseus cluster and concluded no significant 
velocity gradient. In addition, Ota et al. (2007) found that the Suzaku results of the Centaurus cluster 
are difficult to reconcile with claims in Dupke & Bregman (2001b) and Dupke & Bregman (2006). In 
short, previous results by Dupke and Bregman suggest bulk motions in some clusters but with large 
uncertainties. 

Currently the Suzaku XIS (Koyama et al. 2007a) would be the best X-ray spectrometer for 
the bulk motion search, because of its good sensitivity and calibration (Ozawa et al. 2009). In fact, 
Suzaku XIS data were already used for this search in representative clusters. Tight upper limits on 
velocity variations are reported from the Centaurus cluster (1400 km s~^; 90% confidence; Ota et al. 
2007) and A 2319 (2000 km s'^ 90% confidence; Sugawara et al. 2009) among others. 

In order to improve the accuracy of the velocity determination and to search for gas bulk 
motions we analyzed Suzaku XIS spectra of the Abell 2256 cluster of galaxies (A 2256, redshift of 
0.058). This X-ray bright cluster is one of the first systems showing substructures not only in the 
X-ray surface brightness but also in the galaxy velocity distribution (Briel et al. 1991). In the cluster 
central region, there are two systems separating by 3 .'5 in the sky. Motivated by this double-peaked 
structure in their T^OS'Ar image, Briel et al. (1991) integrated the velocity distribution of galaxies from 
Fabricant et al. (1989) over the cluster, fitted it to two Gaussians, and found the two separated peaks in 

' The error confidence range is not explicitly mentioned in the reference. 
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the velocity distribution. The two structures are separated by ~ 2000 km s~ in radial velocity peaks of 
member galaxies, as given by table 1 . Berrington et al. (2002) added new velocity data to the Fabricant 
et al. (1989) sample, used 277 member galaxies in total, and confirmed the two systems along with 
an additional third component (table 1). This unique finding motivated subsequent observations in 
multiple wavelengths. For example, radio observations revealed a diffuse halo, relics, and tailed radio 
emission from member galaxies (e.g. Rottgering et al. 1994). The Chandra observation by Sun et al. 
(2002) revealed detailed gas structures in and around the main and second peaks. Furthermore, there 
are some attempts to reproduce merger history of A 2256 using numerical simulations (e.g. Roettiger 
et al. 1995). Thus, A 2256 is a prototypical and well-studied merging system and hence suitable to 
study the gas dynamics. 



Table 1. Fitting parameters of radial velocity distribution of member galaxies. 



component number Mean velocity 


Dispersion 


in the reference (km s~^) 


(km s-i) 


Brieletal. (1991) 




1 (main) 17880 ±205 


1270 ±127 


2 (sub) 15730 ±158 


350 ±123 


Berrington et al. (2002)* 




1 (53 galaxies; sub) 15700 


550 


2 (196 galaxies; main) 17700 


840 


3 (28 galaxies) 19700 


300 



* Errors on velocities and dispersions are not given in the paper 



We have carefully evaluated instrumental energy scales of the Suzaku XIS, used iron K-shell 
line emission, and found a radial velocity shift of the second gas component with respect to the main 
cluster. The X-ray determined redshifts are consistent with those of galaxy components. This is the 
most robust detection of a gas bulk motion in a cluster. 

Throughout this paper, we assume cosmological parameters as follows; Hq = 70 km 
s~^Mpc~^, fim = 0.3, and = 0.7. At the cluster redshift of 0.058, one arc-minute corresponds 
to 67.4 kpc. We use the 68% (la) confidence level for errors, unless stated otherwise. 

2. Observations 

Suzaku observations of A 2256 were performed on 2006 November 10-13 (PI: K. Hayashida). 
The XIS was in the normal window and the spaced-row charge injection off modes. The observation 
log is shown in table 2. Figure 1 shows an X-ray image of the cluster. Detailed descriptions of the 
Suzaku observatory, the XIS instrument, and the X-ray telescope are found in Mitsuda et al. (2007), 
Koyama et al. (2007a), Serlemitsos et al. (2007), respectively. 
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Fig. 1. Suzaku XIS contour image of A 2256 in the 0.3-10 keV band overlaid on a DSS (Digital Sky Survey by 
the Space Telescope Science Institute) optical image. The three XIS data are co-added. The XIS field of view is 
shown by a square. The contour levels are linear from 10 to 160 counts pixel^^. No vignetting nor background was 
corrected. The calibration-source regions at comers of CCD are excluded. North is up and east is to left. 

To verify the XIS gain calibration, we have also used data from the Perseus cluster observed 
in 2006 with the same XIS modes as those for A 2256 (table 2). These data were already used for the 
XIS calibration (e.g. Ozawa et al. 2009) and scientific analyses (Tamura et al. 2009; Nishino et al. 
2010). 



Table 2. Suzaku observations of A 2256 and the Perseus cluster. 



Target 


Date 


Sequence 


(RA, Dec) 


P.A. * 


Net exposure 






number 


(degree, J2000.0) 


(degree) 


(ks) 


A 2256 


2006/11/10-13 


801061010 


(256.0138,78.7112) 


208 


94.4 


Perseus 


2006/2/1 


800010010 


(49.9436,41.5175) 


260 


43.7 


Perseus 


2006/8/29 


101012010 


(49.9554,41.5039) 


66 


46.6 



* The nominal position angle (north to DETY axis). 



3. Analysis and Results 

3.1. Data Reduction 

We used version 2.1 processing data along with the HEASOFT version 6.9. In addition to 
the XIS standard event selection criteria, we have screened data by adding the following conditions; 
geomagnetic cut-off-rigidity > 6 GV, the elevation angle above the sunlit earth >20° and the dark 
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earth >5°. We used the latest calibration file as of July 2010. Using these files, we corrected the XIS 
energy scale. The data from three CCDs (XIS 0, XIS 1, and XIS 3) are used. 

We examined the light curve excluding the central bright region events (R < 6') for stable- 
background periods. There was no flaring event in the data. The instrumental (non-X-ray) background 
was estimated using the night earth observation database and the software xisnxbgen (Tawa et al. 
2008). 

We prepared X-ray telescope and CCD response functions for each spectrum using software 
xissimarfgen (Ishisaki et al. 2007) and xisrmfgen (Ishisaki et al. 2007), respectively. The 
energy bin size is 2 eV. 

To describe the thermal emission from a coUisional ionization equilibrium (CIE) plasma, we 
use the APEC model (Smith & Brickhouse 2001) with the solar metal abundances taken from Anders 
&Grevesse (1989). 

3.2. Energy Scale Calibration 

The main purpose of this paper is to measure Doppler shifts of K-shell iron lines in X-ray. In 
these analyses, the energy-scale calibration is crucial. In § 3.2.1, we summarise the calibration status. 
In § 3.2.2, we attempt to confirm the calibration using calibration-source data collected during the 
A 2256 observation. The primal goal here is to find differences of line energies observed in different 
positions within the field of view. Therefore the positional dependence of the energy scale is most 
important, which is evaluated in § 3.2.3. Here we focus on the data obtained in spaced-row charge 
injection off mode and around the K-shell iron lines. Considering all available information given here 
as summarized in table 3, we assume that the systematic uncertainty of the energy scale around the 
iron lines is most likely 0.1% and 0.2% at most over the central 14'. 7 x 14'. 7 region or among the 
three CCDs. 
3.2.1. Reported Status 

Koyama et al. (2007b) estimated the systematic uncertainty of the absolute energy in the iron 
band to be within +0.1, —0.05%, based on the observed lines from the Galactic center along with Mn 
Ka and K/3 lines (at 5895 eV and 6490 eV respectively) from the built-in calibration source (^^Fe). 
Independently, Ota et al. (2007) investigated the XIS data of two bright and extended sources (the 
A 1060 and Perseus clusters) and evaluated the positional energy scale calibration in detail. They 
estimated the systematic error of the spatial gain non-uniformity to be ±0.13%. Furthermore, Ozawa 
et al. (2009) examined systematically the XIS data obtained from the start of operation in July 2005 
until December 2006. They reported that position dependence of the energy scale are well corrected 
for the charge-transfer inefficiency and that the time-averaged uncertainty of the absolute energy 
is ± 0.1%. In addition, the gradual change of the energy resolution is also calibrated; the typical 
uncertainty of resolutions is 10-20 eV in full width half maximum. 
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Table 3. Accuracy of the XIS energy scale calibration and the observed velocity shift 



Sectionf 


Line used 


A* (%) 


Reference and remarks 


3.2.1 


GCDX/Fe-Kt 


+0.1,-0.05 


Koyama et al. (2007b); uncertainty of the absolute energy. 


3.2.1 


Perseus/Fe-K 


±0.13 


Ota et al. (2007); spatial non-uniformity (at 68% confidence). 


3.2.1 


MnlKa 


±0.1 


Ozawa et al. (2009); time-averaged uncertainty of the absolute energy. 


3.2.2 


MnlKa 


±0.1 


This paper; standard deviation among CCD segments in the A 2256 data. 


3.2.3 


Perseus/Fe-K 


(-0.09-0.06) 


This paper; variations among the two regions. 


3.4 


A 2256/Fe-K 


0.5±0.1 


This paper; velocity shift in A 2256, after correcting the inter-CCD gain 



and the statistical error 

T Section in this paper. 

* possible uncertainty, velocity shift, or its error all in percentage of the Hne energy. At the iron K hne, 0.1% corresponds to 7 eV for 
the energy shift or 300 km s~^ for the radial velocity, 
t From the Galactic center dift'use emission. 



3.2.2. Absolute Scale 

We extracted spectra of calibration sources which illuminate two comers of each CCD (seg- 
ment A and D). These spectra in the energy range of 5.3-7.0 keV are fitted with two Gaussian lines 
for the Mn Ka and K(3 along with a bremsstrahlung continuum component. Here we fixed the energy 
ratio between the two lines to the expected one. Thus obtained energy centroids of the Mn Ka line 
from the two comers of three CCDs give an average of 5904 eV (as compared with the expected 
value of 5895 eV) and a standard deviation (scatter among the six centroids) of 6 eV. The statistical 
errors of the line center were about 1-2 eV. This confirms that the absolute energy scale averaged over 
all CCDs and the relative gain among CCD segments are within ±0.15% and ±0.10%, respectively. 
Simultaneously we found that the data can be fitted with no intrinsic line width for the Gaussian 
components, meaning that the energy resolution is also well calibrated. 

3.2.3. Spatial Variation 

We used the XIS data of the Perseus cluster, which would provide the highest-quality XIS line 
spectra over the whole CCD field of view among all the Suzaku observations. Here we assume no 
line shift intrinsic to the cluster within the observed region. There were two exposures of the Perseus 
in the normal window and the spaced-row charge injection off modes (table 2) in periods close to the 
A 2256 observation. Note that Ota et al. (2007) used the same data, but with early calibration (i.e. 
version 0.7 data). Here we re-examine the accuracy with latest and improved calibration. 

Firstly, following Ota et al. (2007), we divided the XIS field of view into 8x8 cells of size 
2'.1 X 2M. Each spectra in the 6.2-7.0 keV band is fitted with two Gaussian lines for He-like Ka 
(~ 6700 eV) and H-like Ka (~ 6966 eV) and a bremsstrahlung continuum model. Here we fixed the 
ratio of the central energies between the two lines to a model value of 1 .040 and let the energy of the 
first line be a fitting parameter. Because of the low statistics data in CCD peripheral regions we focus 
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on the central 7x7 cells (i.e., 14'. 7 x 14'. 7). The typical statistical errors of the line energy are from 
±4 eV to ±25 eV. Thus derived central energies from 7x7 cells for each CCD are used to derived 
an average and a standard deviation. The average values are 6575 eV, 6575 eV, and 6569 eV, for XIS 
0, XIS 1 and XIS 3, respectively, which are consistent with the cluster redshift of 0.017-0.018. The 
standard deviations among 7x7 cells are 7 eV, 13 eV, and 10 eV for the three XISs, respectively. 
There is no cell having a significant deviation from the average value at more than 2a level. These 
deviations include not only systematics due to the instrumental gain uncertainty but also statistics 
and systematics intrinsic to the Perseus emission. Therefore the energy scale uncertainty should be 
smaller than this range of deviations (0.1-0.2%). 

Secondly, we focus on the CCD regions specific in our analysis below. We have extracted the 
Perseus spectra from the same detector regions (main and sub) as used in the A 2256 analysis. The 
definition of the regions are given in figure 2 and § 3.4. These spectra were fitted with the same model 
as above and used to derive the redshift (from the line centroid). We found no difference in redshift 
between the two regions. The redshift differences (between the two regions from the three sensors) 
range from —9 x 10""^ to +6 x 10""^ with an average of 1.5 x 10~^. Accordingly we estimate the 
possible instrumental gain shift to be within ±0.1% with no systematic bias between the two regions. 

3.3. Energy sorted X-ray Images 

In order to examine the spectral variation over the A 2256 central region, we extracted two 
images in different energy bands including He-like and H-like iron line emission as shown in figure 2. 
There appear at least two emission components which corresponds to the main cluster at east and 
the sub component at west discovered in Briel et al. (1991). We noticed a clear difference in the 
distribution between the two images. In the He-like iron image, the sub component exhibits brightness 
comparable to that of the main component. On the other hand, in the H-like iron image, the sub 
has lower brightness. This clearly indicates that the sub has cooler emission compared with the 
main. Based on these contrasting spatial distributions, we define centers of the two components to 
be (256.°1208, 78.°6431) and (255.°7958, 78.°6611), in equatorial J2000.0 coordinates (RA, Dec), 
respectively, as shown in figure 2. 

3.4. Radial Velocity Shift 

3.4.1. Separate Spectral Fitting 

Our goal is to constrain the velocity shift of the sub component with respect to the main cluster. 
Then, we extracted sets of spectra from the two components and fitted with different redshifts. More 
specifically, for the main component emission, we integrated the data within 4' in radius from the 
main center but excluding a sub component region with a radius of 2' (figure2). For the sub region, 
data within 1'.5 in radius from the sub center are extracted. 

We use the energy range of the 5.5-7.3 keV around the iron line complex. In this energy band 
the cosmic background fluxes are below a few % of the source count in the two regions. Therefore 
we ignore this background contribution. 
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3.3 4.7 5.8 6.6 7.4 8.1 8.8 9.4 10 

Fig. 2. XIS images of A 2256 from the 6. 10-6.48 keV energy band including He-like iron emission (left) and from 
the 6.48-6.80 keV band including H-like iron emission (right) in unit of counts. No vignetting nor background was 
corrected. Images have been smoothed with a Gaussian filter with cr = 17". The spectra are extracted from the two 
regions indicated by circles in yellow (the main region) and a circle in white (the sub). 

We use a CIE component to model the spectra. Free parameters are an iron line redshift, 
a temperature, an iron abundance, and normalization. Models for the three XISs are assumed to 
have different redshifts and normalization to compensate inter-CCD calibration uncertainties. Other 
parameters, temperature and abundance, are fixed to be common among the CCDs. As shown in 
figure 3 and table 4, these models describe the data well (FIT-1). The fits give different redshifts 
between the two regions from all three XISs as shown in figure 4(a). The statistical errors of the 
absolute redshift are less than ±0.002. The redshift differences are 0.0048, 0.0041, and 0.0076 for 
the three XIS. These are equivalent to shifts of 30-50 eV in energy or 1200-2300 km s~^ in radial 
velocity. 

The spectral shape of the two regions are different as shown in figure 3. Therefore, the obtained 
redshift difference may depend on the spectral modeling. To check this possibility, instead of the CIE 
model, we use two Gaussian lines and a bremsstrahlung continuum component (FIT-2). Here we 
assume the He-like iron resonance line centered on 6682 eV and the H-like iron Lya centered on 
6965 eV for the two lines and determine the redshift common to these two lines. The Gaussian 
components are set to have no intrinsic width. As given in table 4, these models give slightly better 
fits to the data in general than the first model . For each XIS, we obtained redshifts and therefore a 
redshift difference as same as those from the first fit within the statistical uncertainty. Therefore the 
redshift depends insignificantly on the spectral modeling. 

Strictly speaking line centroids of the iron line transitions could change depending on the 
emission temperature. In the case of the observed region of A 2256, where the temperature varies 



8 



within 4-8 keV (e.g. Sun et al. 2002), the strong iron line structure is dominated by the He-like 
triplet. Within this temperature range, the emission centroid of the triplet stays within 6682-6684 eV 
based on the APEC model. This possible shift (< 2 eV) is well below the obtained redshift difference 
(30-50 eV) and should not be the main origin of the difference. 
3.4.2. Gain-corrected Spectral Fitting 

The obtained redshifts are systematically different among the three XISs [figure 4(a)]. We 
attempt to correct this inter-CCD gain difference based on the calibration source data. Following 
Fujita et al. (2008), we estimate a gain correction factor, /gain, by dividing the obtained energy of the 
Mn K line from the calibration source by the expected one. In the A 2256 data as given in section 3.2, 
/gain are 1.0018 (XIS 0), 1.0000 (XIS 1), and 1.0016 (XIS 3). This factor, the redshift obtained from 
the fit (2;fit), and the corrected redshift {Zcor) have a relation 

Zcor + 1 

We use this correction and fit the spectra with a single redshift common to the three XISs to the two- 
Gaussian model (FIT-3). These models give slightly poorer fits compared with the previous ones, 
because of a decrease of the degree of freedom (table 4). Nevertheless, the fits are still acceptable. 
The difference in the redshift between the two region is 0.005 ± 0.0008 (or 1500 ± 240 km s^^). In 
figure 4(b) we show the statistical distribution as a function of the redshift. This indicates clearly 
that the data cannot be described by the same redshift for the main and sub regions. 

We found that redshifts determined here by X-ray are consistent with those of member galaxies 
in optical (table 1) as explained below. The X-ray redshift of the main component, z=0.0059 or 17700 
km s~^, is the same as the galaxy redshift within the statistical error. That of the sub component, 
z=0.00540 or 16200 km s^^, is larger than the galaxy value of 15730 km s^^. Yet the difference, 
470 km s^\ is within the combined errors from X-ray statistical (150 km s^^), systematic (300km 
s~^), and optical fitting (160 km s~^). Besides this, the obtained specta of the sub component could 
be contaminated from the emission of the main component due to the projection and the telescope 
point spread function (half power diameter of about 2'). This contamination could make the obtained 
redshift of the sub larger than the true one (or closer to that of the main). 

To visualize the redshift shift in the sub region more directly, we show spectral fittings by 
combining the two front- illuminated XIS (0 and 3) data in figure 5. We fitted the spectra with the 
two Gaussian model as given above. In this figure we compare the fitting results between the best- 
fit model {z = 0.051) and one with the redshift fixed to the main cluster value (z = 0.058). This 
comparison shows that we can reject that the sub region has the same redshift as the main cluster. 

4. Summary and Discussion 

4.1. Summary of the Result 

We found gas bulk motion of the second component in A 2256. The difference in the redshifts 
and hence radial velocities between the main and sub systems is found to be 1500 ±300 (statistical) 

9 



(a) FIT- 1 , main (b) FIT- 1 , sub 




Energy (keV) 

Energy (keV) 

Fig. 3. Cluster spectra along with the best-fit CIE model in histogram (FIT-1). Plots (a) and (b) are from the 
main and sub regions, respectively. The XIS 0, XIS 1, and XIS 3 data are shown by black, red, and green colors, 
respectively. In lower panels fit residuals in terms of the data to model ratio are shown. 



Table 4. Spectral fitting results from the main and sub regions. 



Region 


kT Fe < z > * 


xVd.o.f. 




(keV) (Solar) 






FIT- 1 , CIE component 




main 


7.1±0.3 0.25±0.02 0.0583 


121/139 


sub 


5.0±0.4 0.29±0.03 0.0528 


135/139 




FIT-2, two Gaussian components 




main 


6.2±0.7 - 0.0580 


112/137 


sub 


4.5±0.9 - 0.0528 


131/137 




FIT-3, two Gaussian with the gain correction 


main 


6.2±0.7 - 0.0590±0.0007 


118/137 


sub 


4.6±0.9 - 0.0540±0.0005 


137/137 


* Average 


value over three XIS data tor the Wl'-l and Hr-2 are 


given. 



Statistical errors on these values are typically 0.001. In the FIT-3, a single 
redshift common to the three XISs is assumed. 
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(a) (b) 




1 2 3 0.055 0.06 

XIS redshift 



Fig. 4. (a) The redshifts obtained from the CIE fitting (FIT-1) of the main and sub component regions for each XIS 
(X-axis). The gain-corrected results of FIT-3 are shown at the position of XIS=1.5. The two horizontal dashed 
lines indicate redshifts from the optical data (Briel et al. 1991; table 1). (b) The statistic distributions of the 
gain-corrected Gaussian fitting (FIT-3) of the two regions as a function of the redshift. The horizontal line at 
= 134 indicates the degree of the freedom of the fits. The two vertical dashed lines indicate the redshifts as in 

(a). 

±300 (systematic) km s~^ (§ 3.4). This observed shift corresponds to only 0.5% in the energy, but is 
well beyond the accuracy of the energy scale reported by the instrument team (Koyama et al. 2007b; 
Ozawa et al. 2009) and that by Ota et al. (2007). Focusing on the present analysis of A 2256, we also 
have examined the calibration systematics independently and confirmed the accuracy given above. 
The obtained redshifts and hence the difference between the two X-ray emitting gas components are 
consistent with those of the radial velocity distribution in member galaxies. This consistency uniquely 
strengthens the reliability of our X-ray measurement. 

This is the fist detection of gas bulk motion not only in A 2256 but also from Suzaku , which 
presumably has the best X-ray spectrometer in operation for the velocity measurement. As given in 
§ 1, Dupke and Bregman (2001a, 2001b, 2006) previously claimed detections of bulk motions in the 
Perseus and Centaurus clusters. Compared with these and other attempts, our measurement is more 
accurate and robust. This improvement is not only due to the better sensitivity and calibration of the 
XIS but also due to the well- separated and X-ray bright nature of the structure in A 2256. 
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(a) FIT-3, sub, best-fit (b) FIT-3, sub, z=0.058 




Energy (keV) Energy (keV) 

Fig. 5. The spectra from the sub-component region by adding the two front-illuminated XIS (0 and 3), fitted with 
two Gaussians and a continuum model, (a) The redshift common to the two lines is left free. The model gives the 
best-fit redshift of 0.051 15 and of 40.4 for the degrees of the freedom of 42. (b) The redshift is fixed to the main 
component value (z = 0.058). The model gives of 80.8 for degrees of freedom of 43. 

Radial velocity distributions of cluster galaxies were used to evaluate dynamics of cluster 
mergers. These studies, however, have limited capability. First, largely because of the finite number 
of galaxies and the projection effect, a galaxy sub group within a main cluster is not straightforward to 
identify. Second, the major baryon in the cluster is not galaxies but the gas in most systems. Therefore 
dynamical energy in the hot gas cannot be ignored. Third, galaxies and gas are not necessary to move 
together especially during merging phases, because of the different nature of the two (coUisionless 
galaxies and coUisional gas). In practice, some clusters under the violent merging phase such as 
IE 0657-558 (Clowe et al. 2006) and A 754 (Markevitch et al. 2003) show a spatial separation between 
the galaxy and gas distributions. Therefore, it is important to measure dynamics of galaxies and gas 
simultaneously. The present result is one of the first such attempts. 

4.2. Dynamics of A 2256 

The determined radial velocity, Vr ~ 1500 km s~^, gives an lower limit of a three dimensional 
true velocity, v = Vr/ sin a, where a denotes the angle between the motion and the plane of the sky. 
Given a gas temperature of 5 keV (§ 3.4) and an equivalent sound speed of 1100 km s~^ around the 
sub component, this velocity corresponds to a Mach number M > 1.4. Therefore, at least around the 
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sub component, the gas kinematic pressure (or energy) can be (1.4)^ ~ 2 times larger than the thermal 
one. In this environment, the gas departs from hydrostatic equilibrium. Then, does this motion affect 
the estimation of the mass of the primary cluster ? As argued by Markevitch & Vikhlinin (1997), it 
depends on the physical separation between the two components. In the case of A 2256, the two are 
not likely too closely connected to disturb the hydrostatic condition around the primary, as estimated 
below. However, to weight the total mass within the larger volume including the sub component we 
should consider not only the mass of the sub itself but also the dynamical pressure associated with 
the relative motion. This kind of a departure from hydrostatic equilibrium was predicted generally at 
cluster outer regions in numerical simulations (e.g. Evrard et al. 1996). 

We will compare our measurement with other studies on the ICM dynamics. Markevitch 
et al. (2002) discovered a bow shock in the Chandra image of IE 0657-56 (the bullet cluster) and 
estimated its M to be 3.0 ± 0.4 (4700 km s~^) based on the observed density jump under the shock 
front condition. Using a similar analysis, Markevitch et al. (2005) derived M of 2.1 ± 0.4 (2300 
km s^^) in A 520. Based on a more specific configuration in the "cold front" in A 3367, Vikhlinin 
et al. (2001) estimated its M to be 1.0 ± 0.2 (1400 km s~^). In these cases, velocity directions are 
assumed to be in the plane of the sky. See Markevitch & Vikhlinin (2007) for a detailed discussion and 
other examples. These measurements are unique but require certain configurations of the collision to 
apply to observations. In contrast. X-ray Doppler shift measurements as given in the present paper 
are more direct and commonly applicable to merging systems. This method is sensitive to the motion 
parallel to the line of sight. These two measurements are complementary and could provide a direct 
measurement of three dimensional motion if applied to a merging system simultaneously. 

We measured the gas velocity parallel to the line of sight (f^). How about the velocity in the 
plane of the sky (Vsky) and the true velocity vl The ROSAT X-ray image of A 2256 shows a steeper 
brightness gradient between the main and second components (Briel et al. 1991). Furthermore, from 
the Chandra spectroscopic data Sun et al. (2002) found a hint of a temperature jump across the two 
component. They argued a similarity of this feature to the "cold front" discovered in other clusters. 
The temperature and hence pressure jump in A 2256 (from about 4.5 keV to 8.5 keV) are similar to 
those found in A 3367, in which the jump indicates a gas motion with fsky of about 1400 km s~^. 
Therefore we expect a similar v^ky in A 2256, which is comparable to v^- Further assuming that the 
second component directs to the main cluster center, v is estimated as ^Jvf+v^ ~ \^Vr ~ 2000 km 
s~^. Instead of assuming Vgky, by considering the mean sin a factor, 2/Tr, v becomes 2400 km s~^ 
on average. Based on a simple assumption that the two system started to collide from rest, we can 
estimate the velocity as f ~ (^^^)^^^, where G,M and R are the gravitational constant, the total 
mass of the system, and the separation, respectively. The total mass of A 2256 of 8 x 10^^ Mq 
(Markevitch & Vikhlinin 1997) and an assumed final separation i? of 1 Mpc give v ~ 2800 km s~^, 
which is comparable but larger than the estimated current velocity. Putting v = 2000 — 2400 km 
s~^ to the relation, we obtain a current separation i? of 1.4 — 2 Mpc. The time to the final collision 
can be estimated to be about 0.2-0.4 Gyr [(R — 1) Mpc divided by v]. Here we assume that the 
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system is going to the final collision. This assumption is consistent with a lack of evidence for strong 
disturbances in the X-ray structure, as argued previously for A 2256 in general. 

4.3. Future Prospect 

Our observations provided the first detection of gas bulk motion in a cluster. To understand in 
general the cluster formation which is dominated by non-linear processes, systematic measurements 
in a sample of clusters are required. For example, some clusters such as IE 0657-558 (Clowe et al. 
2006) at violent merging stages show segregation of the gas from the galaxy and possibility from dark 
matter components. In these systems, we expect different situations in the gas and galaxy dynamics 
compared with that found in A 2256. Given capabilities of current X-ray instruments such as the 
Suzaku XIS, A 2256 is presumably an unique target with X-ray flux high and the velocity separation 
clear enough to resolve the structure. Accordingly, the systematic study requires new instruments 
with higher spectral resolutions and enough sensitivities. In fact, this kind of assessment is one of the 
primary goals for planned X-ray instruments such as SXS (Mitsuda et al. 2010) onboard ASTRO-H 
(Takahashi et al. 2010). Using the SXS with an energy resolution better than 7 eV we could measure 
gas bulk motions in a fair number of X-ray bright clusters. Furthermore, we may find line broadening 
originated from the gas turbulence, as a result of mergers, related shocks, or some activities in the 
massive black hole at cluster centers. In addition, the SXS potentially will constrain for the first time 
the line broadening from the thermal motion of ions (and hence the ion temperature). The present 
result proves that A 2256 is one of the prime targets for ASTRO-H . The expected spectra of the two 
components in A 2256 with the SXS are shown in Fig. 1 1 of Takahashi et al. (2010). 

We thank anonymous referee for useful comments. We thank all the Suzaku team member 
for their supports. To analyze the data, we used the ISAS Analysis Servers provided by ISAS/C- 
SODA. KH acknowledges the support by a Grant-in- Aid for Scientific Research, No.2 1659292 from 
the MEXT. 
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